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Abstract 

We study the universal characteristics of the shape of a polymer chain in an envi- 
ronment with correlated structural obstacles, applying the field-theoretical renor- 
malization group approach. Our results qualitatively indicate an increase of the 
asymmetry of the polymer shape in crowded environment comparing with the pure 
solution case. 
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1 Introduction 



In studying the transport properties of polymer fluids, an important role is 
played by the shape characteristics of a single polymer chain configuration. It 
is established [1,2], that the typical polymer chain realization does not possess 
spherical symmetry. A quantity, which characterizes the asymmetrical shape 
of a polymer chain, is the asphericity ratio Ad of the chain configuration. Ad 
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attains a maximal value of one for a completely stretched, rod-like configu- 
ration, and equals zero for the spherical form (see Fig. 1), thus obeying the 
inequality: < A d < 1. In the limit of long chains this quantity appears to 
be universal and depends on space dimension d only: A d > 1/2 at d < 4, 
A d = 1/2 at d > 4. The size measure of a flexible polymer chain is usually 
defined by either the mean-squared end-to-end distance R e or radius of gyra- 
tion Rq- The ratio of these quantities, the so-called size ratio g = (RI)/(Rq) 
again is an universal quantity (g > 6 for d < 4, g — 6 for d > 4). The study 
of these universal quantities characterizing the polymer shape is a subject of 
a great interest [3,4,5,6]. 




(a) (b) 

Fig. 1. Schematic presentation of a polymer configurations with aspericity ratio A d 
value close to one (a) and close to zero (b). 

A current problem in a polymer physics is the influence of structural obstacles 
(impurities) in the environment on the universal properties of macromolecules. 
Such a disordered (crowded) environment can be found, in particular, in a bi- 
ological cell, composed of many different kinds of biochemical species, which 
may occupy large fraction of the total volume [7]. In the language of lattice 
models of polymers, the crowded environment with structural obstacles can be 
considered as a disordered lattice, where a given fraction of randomly chosen 
sites are to be avoided by the polymer chain. It has been proven both analyti- 
cally [8] and numerically [9] , that weak uncorrelated disorder that corresponds 
to the point-like randomly distributed obstacles with concentration far from 
the percolation threshold does not change the universality class of polymers. 
In the present study, we address a model, where the structural obstacles of 
the environment are spatially correlated on a mesoscopic scale [10]. Following 
Ref. [11], this case can be described by assuming the defects to be correlated 
on large distances r according to a power law with a pair correlation function 
h(r) ~ r~ a . For a < d such a correlation function describes defects extended 
in space (see [11,12] for further details). The impact of long-range-correlated 
disorder on the scaling of single polymer chains has been analyzed in previ- 
ous works [12] by means of the field-theoretical renormalization group (RG) 
approach. The intriguing question of the how the shape characteristics of a 
flexible chain are influenced by presence of such a disordered medium remains, 
however, still unresolved, and is the subject of the present study. 

The layout of the paper is as follows. In the next section, we describe the 
model of the flexible polymer chain in a crowded environment, develop its 
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field-theoretical formulation and give a brief description of the field-theoretical 
RG approach. The results of its application to the present model and finally 
some conclusions and an outlook are given in the following sections. 



2 Field-theoretical description of the model 



Let R n = {x^, . . . , x%} be the position vector of the nth monomer of the poly- 
mer chain (n = 1, . . . , N). The shape of a specific conformation of the chain 
can be characterized [2] in terms of the gyration tensor Q with components: 

1 N 

Qa = w E « - <J« - O, i,j = l,...,d. (1) 

ZJV n,m=l 



For the averaged radius of gyration Rc one thus has: 
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( R o) = ^{ E \Rn-R m \ 2 ) = (T,Q*) = ( 2 ) 



\n,m=l I \i=l 



Here and below, (...) denotes averaging over all configurations of the polymer 
chain. The spread in eigenvalues \ of the gyration tensor measures the asym- 
metry of a given configuration. In Ref. [5] it was proposed to characterize the 
shape of macromolecule by ratios of rotationally invariant polynomials in the 
components of Q. Let us define the asphericity ratio Ad as the quotient of two 
averaged quantities [5]: 

i f ((A.-A) 2 ) d (TrQ 2 ) 



d{d-l)U ((X) 2 ) d-l((TrQy 



Here, A = TrQ/d is the mean eigenvalue of gyration tensor for a given chain 
configuration and Q = Q — AI with the unity matrix I. Ad equals zero when all 
configurations are spherical (eigenvalues of each given configuration are equal) 
and takes the maximal value of one when all configurations are rod-like (all 
the eigenvalues equal zero except one) . Besides the ratio of averages one may 
also calculate the so-called mean asphericity, given by an average of the ratio 
[13]. While ref. [13] has shown some slight quantitative deviations between the 
two approaches it confirmed that the qualitative dependence with respect to 
other parameters, e.g. on the architecture of the molecule is the same in each 
case. 

Passing from the discrete model to a continuous polymer chain description 
[14] allows to derive a field-theoretical formulation of the asymptotic shape 
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characteristics of a single polymer chain in an environment with correlated 
structural obstacles. Applying the replica method in order to average the free 
energy over different configurations of the randomly distributed and fixed 
obstacles leads to an m-component field theory with a Lagrangian L Dis [12]: 



LD- ls = lt J dM(Aol4(x)| 2 + |v4(^)l 2 ) + 5(0l(^)) 2 ] 

a=l 

n „ 

£ J d d xd d yh(\x-y\)$l(x)%(y). (4) 



a,B=l- 



Here each replica (f a is an m-component vector field <p a = . . . ,0™), /to 
and Uq are bare mass and coupling, the coupling of the replicas are described 
by the correlation function h(r) ~ r~ a and both the polymer (m — > 0) and 
the replica (n — > 0) limits are implied. For small k, the Fourier-transform h{k) 
of h(r) behaves as: h(k) ~ v + Wo\k\ a ~ d . Taking this into account, rewriting 
Eq. (4) in the momentum space, and recalling special symmetry peculiarities 
of (4) that appear for m, n — > [12], a theory with two bare couplings Uq, Wq 
results. Note that for a > d the w -term is irrelevant in the RG sense and one 
restores the pure case (absence of structural defects). 

As shown in Ref. [5], computing the shape parameters of long polymer chain 
can be reduced to computing the critical amplitudes and exponents of the 
corresponding field-theoretical model. In order to extract the scaling behavior 
of the model (4) we apply the RG method [15] to get the Green's functions 
G^' N ^ renormalized at non-zero mass and zero external momenta. The change 



of couplings Uq, Wq under renormalization defines a flow in t 
space, governed by corresponding /3-functions: (3 u (u,w^ 



d\n£ 



re parametric 

, (3 w (u,w) = 



dw 



, where / is the rescaling factor, and 

o 



stands for evaluation at fixed bare 

o 



parameters. The fixed points (FPs) u*,w* of the RG transformation are given 
by common zero of /3-functions. The stable FP corresponds to the critical 
point of the system. 

Following Ref. [5] , the averaged moments of the gyration tensor Q needed to 
determine the polymer shape characteristics (2) and (3) can be expressed in 
terms of the connected Green's function. In particular: 



(Qij) = ~\ ( 



DT\ 



2;/ 



{QijQ 



2X J 
1 



r( 7 ) 



G 



r( 7 + 2z/ + 2) Gg } (0,0,{A*}) 



klj 



DT\ 4u 



G 



ij\kl 



A\2XJ r(7 + 4z/ + 4) G g)( ,0,{A*})' 



(5) 
(6) 



4 



where D,T,X are non-universal quantities which will drop out in the final 
expressions, T is the Euler gamma-function, v and 7 are the critical exponents 
of model (4), and the following notations are used: 



G W = f^^J^^r4 2 ' 4) (0,0;q 1 ,-q 1 ,q 2 ,-q 2 ;{A*})') ^ (8) 



dq\ dq{ dq\ dq l 2 



Here Gg' 2) (0, 0; q, -q; {A*}) and Cg' 4) (0, 0; qi, -qi, q 2 , -q 2 ; {A*}) are the fixed 
point values of the renormalized two-point connected Green's functions with 
two and four 2 -insertions, the symbol q indicates that the corresponding 

expressions are to be taken at all external momenta {q} equal to zero. 

The isotropy of the model implies, in particular, that (Tr Q) = d(Q xx ), so 
that: 

(R G )=d(Q xx ). (9) 



For the mean-squared end-to-end distance (R 2 e ) one has [5]: 

(R 2, _ _ (DT\ 2 » M (V^g^-MA*}))!^ 

{ e) \2Xj T{ 1 + 2v) 4 2) (0,0,{A*}) ' { } 



where V^. means differentiation over components of external moment k. The 
not-universal quantities cancel when the ratio g = (RI)/(Rq) is considered: 

_ 2r( 7 + 2u + 2) (VjggCk, -k, {A*})) | k=0 
9 r( 7 + 2i/) (v^(0,0;q,- q ;{A*}))| q . o ( > 



and thus the size ratio is universal quantity. The asphericity ratio A d can be 
expressed in terms of the averaged moments of gyration tensor (5), (6) as 
follows: 

i _ (Qxx) + d(Q xy } — (Q XX Qyy) 

d ~ (QL) + d(d-l)(Q xx Q yy ) ■ [ ) 



One can again easily convince oneself, that all non-universal quantities in Eqs. 
(5), (6) cancel when calculating (12), and Ad is a universal quantity. 
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3 Results 



To obtain the qualitative characteristics of the critical behavior of the model, 
we exploit a double expansion in the parameters e = 4 — d and 5 = 4 — a, 
assuming them to be of the same order of magnitude. Within this approach 
it was shown [12], that a single polymer chain in a solvent in an environment 
with long-range-correlated structural obstacles belongs to a universality class 
different from the case of a pure solvent. In the field-theoretical renormaliza- 
tion group description, this is reflected by the appearance of a new non-trivial 
stable long-range-correlated (LR) fixed point besides the usual (pure) one. 
The coordinates of these FPs and their regions of stability read [12]: 



pure FP : u* = f , w* = 



LR FP : u* 



4 

3S 2 
2(«-e)' 



W 



at 5 < e/2, 
ffEf ate/2<5<e. 



(13) 



To estimate the size ratio (11) for the case of a polymer in long-range-correlated 

(2 2) 

disorder, we calculate the function G R ' (0, 0; q, — q; u, w) with two insertions 
2 (q), 4> 2 (— q). In the first order of the e = 4 — d, 5 = A — a-expansion we find: 



2 4 1 2 
Gg' 2) (0,0;q, -q;u,w) = - -^-—[uh-wJ^ - - [uI 2 -wJ 2 } 

q z + 1 6 q z + 1 6 

+ t^TT [uI °~ wJo] - (14) 

Here Jj, Jj are given by the following one- loop integrals: 



_ f dp _ f dp\p\ a d 

Io -J p+ip Jo " J WTTy (15) 

r f dp r dp\p\ a - d 

1 J (p 2 + l)((p + g) 2 + l)' 1 J (p 2 + l)((p + g) 2 + l)' 1 J 

r = /" d P j = [ dp\p\ a ' d 

2 J (p 2 + l) 2 ((p + g) 2 + l)' 2 7 (p 2 + l) 2 ((p + g) 2 + l)' 1 ; 

To evaluate g according to Eq. (11), we perform the derivation of G R ' with 
respect to the components of the vector q and expand the loop integrals in e 
and 5. Further inserting the FP values (13) we finally receive: 

P urc = 6 + jL S<e/2, 

16 ' (18) 

g ^ = 6 + I , e/2 < 5 < e. 
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Let us qualitatively estimate the change in the size ratio g, caused by presence 
of the structural obstacles, in three dimensions. Substituting directly e — 1 
into the first line of (18), we have for the polymer chain in a pure solvent: 
gpurc ^ g Qg us recall, that the influence of the long-range-correlated 
disorder is relevant for a < d (5 > e) (see e.g. explanation after Eq. (4)). Esti- 
mates of g LR can be evaluated by direct substitution of continuously changing 
parameter S into the second line of Eq. (18). One concludes, that increasing 
the parameter 5 (which corresponds to an increase of strength of disorder) 
leads to corresponding increase of the ratio of the end-to-end to the gyration 
radii g. 

To compute the averaged asphericity ratio using (12), we calculate the function 
G^' 4) (0,0; qi, -qi,q 2 , -q 2 ;u,w) with four insertions 2 (qi), 2 (-qi), 2 (q2), 
2 (— q 2 ). The resulting expansion is too cumbersome to be presented here 
and will be given elsewhere [16]. Performing the derivatives with respect to 
the components of the vectors qi, q 2 we find: 



G xx = 576 + —— (it - w), 
15 

436 

G xx \yy = 320 + — (it - W), 

308 

G X y\ X y = 128 + —(U - W). 



(19) 



Inserting the FP values (13) into Eqs. (19) and recalling the definitions (5) 
and (6) the result is: 



pure 



A d 
t L R 



2 ~ 512 c > 



5 < e/2, 



l + i^ + IH e/2<5<e. 



(20) 



Again, let us qualitatively estimate the change in A d caused by the presence 
of structural obstacles in three dimensions. Substituting directly e — 1 into 
the first line of (20), we have for the pure case: A d ~ 0.53. Estimates 

A LR, 

of A d can be obtained by direct substitution of the continuously changing 
parameter 5 into the second line of Eq. (20). Increase of the strength of disorder 
correlations results in increase of the asphericity ratio of polymers in disorder. 
This phenomenon can be intuitively understood if one recalls an impact of the 
long-range-correlated disorder on the mean end-to-end distance exponent v. 
Indeed, it has been shown in [12], that such a disorder leads to an increase of u, 
and subsequently, to further swelling of the polymer chain. Extended obstacles 
disfavor return trajectories and as a result the polymer chain becomes more 
elongated. This elongation then leads to an increase of the asphericity ratio 
as predicted by Eq. (20). 
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4 Conclusions and outlook 



The universal characteristics of the average shape of a polymer coil configura- 
tions in a porous (crowded) environment with structural obstacles have been 
analyzed considering the special case, when defects are correlated at large dis- 
tances r according to the power law: h(r) ~ r~ a . Applying the field-theoretical 
RG approach, we estimate the size ratio g = (R^) / (Rq) an d averaged aspheric- 
ity ratio Ad up to the first order of a double e = 4 — d, 5 = 4 — a expansion. 
We have revealed, that the presence of long-range-correlated disorder leads 
to an increase of both g and A^ as compared to their values for a polymer 
chain in a pure solution. Moreover, the asphericity ratio A^ value was found 
to be closer to the maximal value of one in presence of correlated obstacles. 
Thus, we conclude, that the presence of structural obstacles in an environ- 
ment makes the polymer coil configurations to be less spherical. Let us note 
that the present results, obtained in the first order of an e, 5-expansion should 
serve as a qualitative estimate rather than an accurate numerical evaluation. 
The next step in our analysis will be to obtain the higher order expansions 
for the quantities of interest and to perform computer simulations in order to 
confirm these theoretical results by numerical values. Further details of our 
calculations as well as a generalization of the presented results to the case of 
polymers with complex topology will be given elsewhere [16]. 
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